Coherence Time of a Solid-State Nuclear Qubit 
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We report NMR experiments using high-power, RF decouphng techniques to show that a^^Si 
nuclear spin qubit in a sohd silicon crystal at room temperature can preserve quantum phase for 
10^ precessional periods. The coherence times we report are longer than for any other observed 
solid-state qubit by more than four orders of magnitude. In high quality crystals, these times 
are limited by residual dipolar couplings and can be further improved by isotopic depletion. In 
defect-heavy samples, we provide evidence for decoherence limited by 1// noise. These results 
provide insight toward proposals for solid-state nuclear-spin-based quantum memories and quantum 
computers based on silicon. 
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Quantum information processing devices outperform 
their classical counterparts by preserving and exploit- 
ing the correlated phases of their constituent quantum 
oscillators, which are usuaUy two-state systems called 
"qubits." An increasing number of theoretical proposals 
have shown that such devices allow secure long-distance 
communication and improved computational power 
Solid-state implementations of these devices are favored 
for reasons of both scalability and integration with exist- 
ing hardware, although previous experiments have shown 
limited coherence times for solid-state qubits. The devel- 
opment of quantum error correcting codes 'ij and fault 
tolerant quantum computation 3] showed that large- 
scale quantum algorithms are still theoretically possible 
in the presence of decoherence. However, the coher- 
ence time must be dauntingly long: approximately 10''' 
times the duration of a single quantum gate, and prob- 
ably longer depending on the quantum computer archi- 
tecture jlj. The question of whether a scalable imple- 
mentation can surpass this coherence threshold is not 
only important for the technological future of quantum 
computation, but also for fundamental understanding of 
the border between microscopic quantum behavior and 
macroscopic classical behavior. 

Experimentally observed coherence times (T2) for var- 
ious qubit implementations are shown in Table ^ The 
atomic systems shown — trapped ions and molecular 
nuclei in liquid solution — have already been employed 
for small quantum algorithms ^'^^ coincidentally, 

they show very high values of Q, the product of the 
qubit frequency ojo/Stt and 'kT2. Most solid-state qubits 
show smaller values of Q; as we demonstrate in this Let- 
ter, however, the long coherence times we observe in 
solid-state ^^Si nuclei afford them a Q and VIT2 as high 
or higher than atomic systems, indicating this system's 



promise for solid-state quantum computing. 

Many promising qubits and coherence time measure- 
ments are not mentioned in TableQleither because reliable 
experimental data are not available, or because the ex- 
isting experimental data are taken under conditions not 
sufficiently similar to the corresponding quantum com- 
puter architecture. For example, electron spins bound to 
phosphorous donors in pure, isotopically depleted silicon 
also show promisingly long coherence times 0]. How- 
ever, quantum computer architectures employing semi- 
conductor electron qubits |15i. il6j require semiconductor 
alloying and metallic nanostructures for wave-function 
engineering. These additions, rather than bulk phenom- 
ena, are likely to dominate decoherence for these qubits. 
Likewise, Kane has proposed an architecture for silicon- 
based quantum computing in which embedded ^"'^P nu- 
clear spins are cooled, controlled, and measured with the 
engineered wave functions of donor electrons • Charge 
fluctuations couple to the nuclei via the hyperfine inter- 
action, leading to a nuclear decoherence source not neces- 
sarily observed in bulk experiments. However, we believe 
that bulk solid-state nuclei such as those studied here can 
be used for quantum computing without strongly per- 
turbing their local environment. Reference lid shows one 
possible architecture for quantum computing in nearly- 
pure, bulk crystalline silicon, motivating our measure- 
ments to explore limits in the coherence time of bulk 
^^Si nuclei. 

The qubit in this study is a thermal ensemble of N 
^^Si nuclei, which are spin-1/2, in a single crystal of sil- 
icon. All other isotopes of silicon are spin-0. If the ^^Si 
nuclei were uncoupled and controlled by identical mag- 
netic fields, the system would have only two accessible 
degenerate energy levels. However, the dipolar coupling 
lifts the degeneracy and takes superpositions of the two 
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TABLE I: Experimental Measurements of Coherence Times T2 for Various Qubits. These times were all measured by spin-echo 
Ramsey spectroscopy or four- wave-mixing, except where noted. Q is the product of the fundamental qubit frequency cjq /Stt 
and nT2. The flTa column shows the experimental product of the Rabi frequency and the coherence time; this provides a 
more realistic measure than Q of the available number of sequential single-qubit gates. Technical improvements can increase 
f2T2 as far as Q in some architectures, while for others Q. must be limited in order to maintain selective qubit control. The 
JT2 column shows the product of the coherence time and the measured or expected qubit-qubit coupling speed; roughly the 
available number of sequential two-qubit gates. 
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A recent measurement employing continuous measurement of Rabi oscillations indicates a T2 of 2.5 /xs in this system 0. 
''This result used the decay of low-visibility Rabi oscillations. 
'^As shown in this Letter. 



qubit states into much larger, entangled superpositions 
of all the 2^ available states. These dynamics are nearly 
independent of temperature. The qubit information is 
not completely lost by this dipolar scrambling, since in- 
formation only leaks from the 2^ nuclear states to the 
lattice on a timescale much slower than the dipolar dy- 
namics. The quantum information can be "refocussed" 
by decoupling pulse sequences such as WAHUHA 
and its offspring, which periodically reverse the dipolar 
evolution, extending T2 by several orders of magnitude. 

The nuclear dynamics during such sequences are usu- 
ally described using average Hamiltonian theory (AHT), 
in which the effects of the periodic RF field on the system 
Hamiltonian are expanded in powers of ic, the sequence 
period |l9j . The system Hamiltonian for a nuclear dipo- 
lar solid in the rotating frame may be written as 

n = -Y. nw,i^ - J2 [I, • Ik ~ 3/|/|] , (1) 

where \j is the spin vector and LOj the shift of the Lar- 
mor frequency from the RF carrier for the j**^ nucleus. 
The dipolar coupling coefficients are Djk = ?i^7^(l — 
Scos^ 9 jk)/2r J)., where rjkCOsOjk = {vj — r^) • z. The 
values of ujj and Djk are different for each nucleus due to 
the inhomogeneous magnetic field and the random dis- 
tribution of the ^^Si nuclei in the crystal lattice. The 
sequence we employ in this study consists of 16 prop- 
erly phased and separated 7r/2 pulses, forming MREV- 
16, a variant of the MREV-8 sequence (20(|. With perfect 
pulses, the MREV-8 sequence results in the zero-order 
average Hamiltonian H^°') = -{l/S)T,j ^^ji^j ^ Ij)^ 
where the sign of the x-component of the effective field 
depends on the helicity of the sequence. The MREV- 
16 sequence, shown in Fig. ^ concatenates each helicity 
leading to n^°^ = -Ej(^j/3)/^ Although MREV- 



16 has reduced spectroscopic resolution over MREV-8, 
maintaining the effective field in the z-direction allows 
easier nuclear control. The inhomogeneous offsets de- 
scribed by H^°^ cause static dephasing, which we peri- 
odically refocus by applying 7r-pulses every 120 cycles 
of MREV-16. We employ the Carr-Purcell-Meiboom- 
Gill (CPMG) phase convention to correct for 7r-pulse 
errors 21], as shown in Fig. ^ Such inserted pulses 
would likewise be employed in a NMR quantum com- 
puter for decou plin g and recoupling of multiple dipolar- 
coupled qubits |22| ■ We note that the CPMG convention 
changes the sequence according to the initial phase of the 
nuclear qubit, which is incompatible with the memory of 
unknown or entangled quantum states. However, more 
complex NMR techniques such as composite pulses 
are known to allow phase-independent pulse correction. 
The MREV-16 sequence and the observed T2 times are 
otherwise independent of the initial nuclear phase. 

The experiments were performed using a 7 T su- 
perconducting solenoid NMR magnet (Oxford) and a 
commercial NMR spectrometer (Tecmag) operating at 
59.575 MHz. The probe was homemade, employing vari- 
able vacuum capacitors (Jennings) and high-power ce- 
ramic capacitors (HEC) to accommodate the pulse pow- 
ers of approximately 700 W while still allowing flexible 
tuning for frequency and impedance matching as well as 
phase transient symmetrization. The RF was detuned 
about 120 Hz from the center of the nuclear resonance 
frequency. Pulsed spin locking E^ ] was observed when 
TT-pulses were applied every 5-10 cycles and in spin-echo 
experiments without decoupling, but this effect dimin- 
ished as the TT-pulse frequency was reduced. Each echo 
time-series was taken in a single measurement after ther- 
mal equilibration for half to five times Ti. 

As shown in Fig. [H this combined "CPMG-MREV- 
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FIG. 1; Schematic of the CPMG-MREV-16X 120 pulse sequence. The echoes shown in the upper left and expanded in the 
upper right are data from an isotopically natural single crystal of silicon. These are obtained by first exciting the sample with 
a 7r/2-pulse of arbitrary phase 0, decoupling with the MREV-16 sequence shown in detail on the bottom line, and refocussing 
with TT-pulses of phase = <^ + 7r/2 every 120 cycles of MREV-16. The magnetization is sampled once per MREV-16 cycle in 
the windows marked with an S. 



16x120" sequence allows the observation of hundreds of 
spin-echoes. In a cylindrical sample of single-crystal sil- 
icon with isotopic content 96.9% ^^Si, the T2 due to un- 
refocussed dipolar-couplings is 450 ^s for the [0011 orien- 
tation, as reported previously for this sample [25j. The 
CPMG-MREV-16xl20 sequence extends the coherence 
time in this sample to nearly 2 seconds. This experi- 
ment used a 6 cm-long, 1.5 cm coil wound using 2 mm- 
diameter wire with variable pitch to improve RF homo- 
geneity. The 7r/2-pulse duration was 9.6 /iS for this coil. 
In isotopically natural silicon, the coherence time is even 
longer due to the scarcity of ^^Si in the lattice, as ^^Si 
is naturally only 4.7% abundant. We examined a high- 
quality single-crystal sample (Marketech) with less than 
5 X 10^"^ cm~^ n-type impurities. It was cut into a sphere 
of diameter 1.5 cm and fit tightly in a constant-pitch coil 
approximately 6 cm long. The 7r/2-pulse duration was 
15 /xs for this coil. As shown in Fig.^ the spin-echoes in 
this sample last for as long as a minute, showing a T2 of 
25.0 ±0.2 sec. 

MREV-16 is expected to reduce the effects of the dipo- 
lar coupling to second order in tc by AHT, even in the 
presence of pulse errors Evidence for this order of 
reduction comes from the observation that T2 increases 
quadratically as tc is decreased (Fig. 12)). The spin-echo 
data was analyzed by Fourier-transforming each echo and 
following the decay of the detuned side-peak, integrated 
between half-maxima. A center peak at the RF frequency 
usually appeared in the spectra due in part to ringdown 
effects and in part to nuclear spin-locking to an effec- 
tive field, which tips slightly away from the z-axis in 
the presence of pulse imperfections. The least-squares 
fit to a single exponential was always used to extract the 
value of T2. The decrease in T2 for small tc occurs be- 
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FIG. 2: Coherence time vs. cycle time. The solid line is 
a fit showing the exponent —2.09 ± 0.07 for the isotopically 
enhanced sample (left) and —2.00 ± 0.2 for the isotopically 
natural sample (right). The insets show the integrated log- 
magnitude of the spin-echoes decaying in time for a few cycle 
times. 



cause MREV-16 begins to fail when the ratio of the pulse- 
width to the time t between pulses approaches unity |l9l| . 
The maximum T2 observed (25 sec) is not a fundamen- 
tal result; it may be exceeded using higher-order pulse 
sequences or isotopic depletion. 

As the strength of the dipolar coupling is further de- 
creased by isotopic depletion, the dipolar coupling con- 
stants Djk become much smaller than the frequency 
offsets ujj. The dominant second-order dipolar aver- 
age Hamiltonian term leading to decoherence is then 
the dipolar/offset cross term, which scales as Ti.^^^ cx 
tclZ^jfellw^l. For isotopic percentage p less than about 
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FIG. 3: Echo decay curves for pure poly crystalline silicon 
of natural isotopic abundance. The solid line is a fit to a 
double exponential with time constants 1.6 and 9.8 seconds. 
No significant variation in any timescale is observed as tc is 
changed. 



required to effectively decouple the ubiquitous nuclear 
spins in this system may be inconsistent with millikelvin 
operation, even if small, high-Q coils and low power, win- 
dowless sequences ji^ are employed. For this reason, we 
believe isotopically depleted silicon to be a more promis- 
ing material for nuclear quantum memory, assuming that 
efficient methods for transferring quantum information to 
and from its well-isolated nuclei can be found. 

This study was sponsored by the DARPA-QuIST pro- 
gram. T.D.L. was supported by the Fannie and John 
Hertz Foundation. We acknowledge N. Khaneja for use- 
ful discussions. 



10%, we would expect T2 to be proportional to the dipo- 
lar coupling constants, which vary as the inverse cube of 
the distance between ^^Si isotopes. Correspondingly, we 
expect Tj"^ approaching T^^ as p — > 0, indicating 
that the threshold for fault-tolerant quantum computing 
could be obtained with reasonably depleted silicon. 

We did observe another decoherence mechanism that 
could appear before the Ti limit in real silicon crystals. 
In a polycrystalline sample with shallow impurity content 
similar to the single-crystal sample, we found the T2 to 
be limited to a constant value of about 8 seconds, inde- 
pendent of tc (Fig. O . We also examined a p-type single 
crystal sample with roughly 1-1.5% ^^Si. This sample 
was known to contain a substantial number of disloca- 
tions in addition to 5 x 10"'^^ cm~'^ aluminum impurities. 
Although the data were noisy due to the reduced signal, 
T2 was observed to be no longer than this same timescale 
of 8 seconds, even after 10-fold averaging. 

We believe this reduced T2 in polycrystalline silicon is 
due to the same low-frequency noise source that leads to 
1/ f noise in silicon wafers |2^. This noise is attributed to 
charge traps at lattice defects and other deep impurities, 
which lead to slow fluctuations of the diamagnetic shield- 
ing seen at nearby nuclei. Free carriers and fixed dipolar 
paramagnetic impurities are unlikely to be responsible 
for this T2, since these lead to isotropic magnetic noise 
with correlation times much shorter than the Larmor pe- 
riod. It follows that T2 due to these sources would be on 
the same order of magnitude as Ti, but we measured Ti 
to be nearly 5 hours, consistent with earlier studies [27l |. 
The elimination of 1// noise from oxides and interfaces 
poses a critical fabrication challenge in quantum comput- 
ing designs based on semiconductor impurities 16j 17j 
and Josephson junctions 0,0,^3' t>ut this noise is ex- 
pected to be very small in high-quality bulk single-crystal 
silicon at low temperature. 

Decoupling pulse sequences have been proposed for 
nuclear memory in high mobility GaAs/AlGaAs het- 
erostructures j23|. We caution that the large RF power 
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